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Abstract 

We show that spatial entanglement of two twin images obtained by parametric down-conversion 
is complete, i.e. concerns both amplitude and phase. This is realised through a homodyne de- 
tection of these images which allows for measurement of the field quadrature components. EPR 
correlations are shown to exist between symmetrical pixels of the two images. The best possible 
correlation is obtained by adjusting the phase of the local oscillator field (LO) in the area of max- 
imal amplification. The results for quadrature components hold unchanged even in absence of any 
input image i.e. for pure parametric fluorescence. In this case they are not related to intensity and 
phase fluctuations. 



1 



I. INTRODUCTION 



Optical systems which display quantum entanglement properties in the spatial domain are 
of great interest for applications, since the amount of information that can be manipulated 
and processed in parallel exploiting quantum correlation effects increases substancially with 
respect to the case of single mode beams. Recently, there has been a rise of interest in the 
utilisation of entangled beam in optical imaging (quantum imaging) [1], 0, |^, |J]. 

In this paper, we consider the field generated through the process of frequency down- 
conversion in a travelling wave optical parametric amplifier (OPA). In [TJ, Q it was demon- 
strated that such a system, when coupled with an appropriate classical imaging device, 
is able to generate two symmetrical amplified copies of an injected input image that are 
strongly correlated one to each other: they indeed display synchronized local intensity fluc- 
tuations at the level of quantum noise, and for this reason they may be referred to as twin 
images. 

Here we present new results, that consolidate and complete the picture, showing that the 
two output images are locally correlated, not only with respect to intensity fluctuations, but 
also to "phase" fluctuations. To carry out a phase-sensitive measurements, we consider a 
homodyne detection scheme that allows us to compare the fluctuations of field quadrature 
components from two corresponding pixels of the two output images. We find in general 
that, for an arbitrary quadrature component characterized by the phase 0^ of the local os- 
cillator, the difference between the fluctuations measured in two symmetrical pixels displays 
exactly the same spectrum as the sum of the fluctuations in the orthogonal quadrature com- 
ponent <pi + 7r/2. The common value can be reduced well below the shot noise level over 
the whole image area, provided that the amplification is large enough and the phase (f>L is 
correctly adjusted. The choice of the phase is crucial to obtain a large level of quantum 
correlations between symmetrical pixels for the quadrature component <pL and anticorrela- 
tion for the quadrature component <pL + tt/2. Thus, the twin images exhibit a complete 
Einstein- Podolsky- Rosen (EPR) || entanglement with respect to continous variables ||. 

Since quantum correlations are shown to exist for any couple of symmetrical pixels over 
the whole area of the output images, we speak of spatial quantum entanglement. The system 
exhibits a spatial realization of the EPR paradox for two orthogonal quadrature components 
of the output field similar to that shown in [|7| for the case of the parametric oscillator below 



2 



threshold. In comparison with the analysis of 0, we consider here also the case in which 
a coherent image is injected into the system. Furthermore, the consideration of an OPA, 
instead of an optical parametric oscillator with spherical mirrors, allows for obtaining here 
completely analytical results. 

In presence of an input image, the mean output field is different from zero, and therefore 
amplitude (i.e. intensity) and phase fluctuations correspond to special cases of quadrature 
fluctuations. Therefore the previous analysis allows for concluding in a rather straightfor- 
ward way that in symmetrical pixels not only quantum intensity fluctuations are strongly 
correlated, but phase fluctuations are anticorrelated in the same amount. 

The paper is divided as follows. After a presentation of the optical image amplification 
scheme in section II, in the third section we study the fluctuation spectrum of the quadra- 
ture components measured in homodyne detection. The fourth section is devoted to the 
discussion of amplitude and phase fluctuations. The final section includes conclusions and 
perspectives. 

II. OPTICAL IMAGE AMPLIFICATION SCHEME 

The experimental procedure to generate a pair of quantum entangled images through the 
process of parametric down-conversion close to the degenerate frequency has been discussed 
in previous papers || |J and can be summarized as follows. The x crystal is enclosed 
between two lenses L and L', as shown in Fig.I. We take the z axis as the main light prop- 
agation direction and indicate with x = (x, y) the point coordinates in a generic transverse 
plane. Not shown in the figure is the coherent pump field that activates the process of down- 
conversion and which we take as an ideal classical monochromatic plane wave of frequency 
u p , propagating inside the crystal along the z direction. 

The crystal slab of width l c , ideally infinite in the transverse directions, is cut for type I 
quasi-collinear phase- matching at the degenerate frequency u p /2. Under these assumptions, 
each elementary down-conversion process corresponds to the splitting of a pump photon of 
frequency uj p into a pair of photons of frequencies oj p /2 + Q and oo p /2 — Q (with Q <C u p ), 
propagating with the same polarization and with opposite transverse wavevectors q and —q. 

We designate by ax(x,t), a 2 (x,t), a 3 (x,t), a^(x,t), the slowly varying envelope operators 
of the down-converted field (with respect to the carrier frequency uj p /2) in the input plane 
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Pi, the entrance plane P 2 of the x crystal slab, its exit plane P 3 , and the output plane P 4 , 
respectively (Fig. 1). We denote by cii(x, fl), a>i(q,Q) {i = 1,..4) their Fourier transforms 
in time and in space-time respectivelly. The purpose of the two lenses is to map the Fourier 
plane (q x , q y ) into the physical plane (x, y). In this manner, if an optical image is injected at 
the degenerate frequency oj p /2 in the object plane Pi, the system amplifies portions of this 
image rather than a band of its g-vectors. Indeed, the input-output transformation which 
describes propagation inside the crystal in the linear regime, assuming that pump depletion 
and losses are negligible, can be written as |J: 

a 3 (q, fl) = u(q, fl)a 2 (<f, fl) + v(q, fl)a\(-q, -fl) , (1) 

The presence of the lenses converts it into a relation between the real-space field operators 
in the plane Pi and P4: 

a 4 (x, fl) = u(x, 0)ai(— x, fl) +v(x,fl)a\(x,—Q) , (2) 

where 

u(x,Q) = -u , v(x,Q) = v ^^pflj , (3) 

/ is the focal length of the two lenses and A is the wavelength of the down-converted field. 
The explicit expressions of the gain coefficients u(q,Q) and v(q,Q) can be found in ||. 
Here we just notice that they depend on the linear gain parameter a and on the dispersion 
properties of the crystal; they are functions of the modulus of q and Q and, for Q = 0, display 
a broad maximum, corresponding to that transverse wavenumber which is phase-matched 
at the degeneracy frequency, for which (within the paraxial approximation) 

q 2 = k] - (k p /2) 2 « k s Ao , (4) 

where k s , k p are the wave numbers of the signal and the pump field at the carrier frequencies 
LVp/2 and u p , respectivelly, and A = 2k s — k p is the collinear phase mismatch parameter 
which is assumed non negative. The width of the plateau around the value (H) is on the 



order of go = yk s /l c , the variation scale of \u\ and \v | in the spatial frequency domain. 

We underline that all the results that follow do not depend on the particular form of the 
gain functions, but rely on the fact that they satisfy the following unitarity conditions 

\u(q,Q)\ 2 -\v(q,n)\ 2 = 1, (5) 
u(q, fl)v(— q, — fl) = u(— q, — fi)f (q, fl), 



which guarantee that the free field commutation rules are preserved: 

'(H{q,n),a\tf,tfj\ = 6tf-?)6(Sl-tf), (6) 
[ai(q,t),aitf,{)] = 0. (i = 1,2,3,4) 

On the other hand, with respect to other systems which exhibit input/output relations of 
the same form (e.g. optical parametric oscillators, see e.g. [|Kj), the large spatial bandwidth 
go of the amplifier makes this travelling-wave scheme a good candidate for high resolution 
image amplification. 

For the scheme of Fig.l, the region in the transverse plane which can be efficiently am- 
plified without distortion has a linear size on the order of 

which represents the width of the plateau of the real-space gain functions @. Such a region 
has either the shape of a disc of area ~ Sq = x\ centered at the origin, or a ring of width 
~ Xq, depending on the possibility to have collinear (A = 0) or non-collinear (A > 0) 
phase-matching at Q = 0, respectivelly. We assume that the input image is a coherent 
stationary field of frequency u p /2 confined in this region of plane Pi (see Fig.l) so that 

(ai(f,fl)> = V2^6(fl)a in (x) . (8) 

As explained in details in @, || , whenever the input image is symmetric with respect to the 
system axis, the device works as a phase-sensitive amplifier (see in this connection also [|TT|j). 
In this case, the phase of the input image must be selected in order to optimize the gain. 
Assuming the input image is duplicated before amplification by means of a classical imaging 
device which allows to obtain a symmetrical field ditribution (i.e. ai n (—x) = cti n (x)), the 
system is able to generate in the output plane two amplified copies that are far better 
correlated in space-time than the originals, meaning by this that they display perfectly 
(in the ideal case) synchronized local intensity fluctuations. It was also demonstrated Q 
that in the limit of high gain, the signal-to-noise ratio as measured from a small portion 
of the input image before duplication is preserved in the corresponding portions of the two 
output images: noiseless amplification is therefore achieved for both output channels taken 



separately (see |12| for an experimental observation of noiseless amplification of images). 

In H an alternative way to generate a pair of quantum correlated images (also called 
twin images) was considered; it consists in injecting a single input image asymetrically, for 



example by confining it to the upper half of plane Pi as shown in Fig.l. This configuration 
does not require a duplication system and presents the further advantage that the gain does 
not depend on the phase of the input field because the systems works as a phase-insensitive 
amplifier. However, the fidelity with which information is transferred is worse than in the 
phase-sensitive case, since the signal-to-noise ratio is deteriorated at least by a factor two in 
the amplification process (a feature common to all phase- insensitive optical amplifiers [|TB||). 

Most of the results presented in this paper do not depend on the particular injection 
scheme, so no assumption are made on the input intensity distribution |a in (x)| 2 . Imperfect 
detection can be modelled in the usual way, by coupling the output field operator a^x, t) with 
an independent operator field a^{x, t) which acts on the vacuum state. The contribution 
describes the noise added by losses in the detection process; thus the effective output field 
measured by a detector of quantum efficiency rj < 1 is 

a D (x, t) = ^/rja 4 (x, t) + ^1- rja N (x, t) (9) 

and the corresponding photon flux density is 

i{x,t) = a} D {x,t)a£>{x,t) . (10) 

As shown in Fig.l, at the exit face of the crystal we insert a pupil of area S p , an element that 
allows for eliminating divergencies which arise in the calculation of the field mean intensity 
and correlation functions when dealing with a system of infinite transverse dimensions H, 



T0|| • It also determines the characteristic resolution area of the device in the detection 
plane, which is Sr = (\f) 2 /S p . This finite size optical element introduces a convolution 
integral with the pupil response function in the r.h.s. of Eq. (|2]) and, as a consequence, the 
points of the input image are spread into diffraction spots of area Sr in the output image. 
However, analytical calculations are performed in the low diffraction limit, assuming that the 
diffraction spot size \/~Sr is much smaller than both xq and the variation scale of the input 
image intensity. Considering a single pixel detector (labelled by index j) that intercepts the 
photons arriving on an area Rj which is large in comparison with Sr, the mean value of the 
measured photocurrent is then |§]: 

= dx {i(x, t)) 

f Tj f f 00 dQ! r 

= r] dx \u(x, 0)a in (-x) + v(x, 0)a* n (f)| 2 + — dx — — \v(x, f2')| 2 . (11) 

JRj Or JRj J-oo 2,71 
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The first integral represents the amplified coherent input field while the second integral is 
the contribution coming from spontaneous parametric down-conversion. The ratio Sq/Sr 
gives an evaluation of the number of details of the input image which can be resolved in the 
detection plane (e.g. with the pixel array of a CCD camera). Moreover, quantum correlation 
effects tends to disappear when Sr — > So, since in this limit the signal and idler photons of 
each down-converted pair can no more be resolved separately, because of the large diffraction 
spread in g-space. Making Sr as small as possible with respect to So is therefore a necessary 
requirement that must be taken into account in experiments. However, this leads to an 
increase of the spontaneous emission contribution which goes at the expense of the visibility 
of the amplified input image. This last circumstance imposes a lower limit on the intensity 



of the input image (see B |8|, [TO] for more details). 



III. CORRELATIONS MEASUREMENT IN A HOMODYNE DETECTION 
SCHEME 

A homodyne detection scheme allows for the measurement of a particular quadrature 
component of the field. It consists in a beam splitter that combines the output field with a 
coherent field of much higher intensity, «l(x), which can be treated as a classical quantity 
and is usually referred to as the local oscillator field (LO). In the balanced version a 50/50 
beam splitter is used, so that the operators associated to the fields coming from the two 
output ports of the beam-splitter , labelled by b and c, are 

a b > c (x , t) = [a 4 (x, t) ± a L {x)\ /V2 (12) 

and the effective fields seen by two identical detectors of quantum efficiency r] in the two 
ports b and c are 

4 c (f , t)=y/fj a»>%x, t) + a% c (x, t) , (13) 

where a^{x,t) describe the noise added in the detection process. When the corresponding 
intensities are electronically substracted, one obtains a direct measure of the quadrature 
component of the output field a 4 selected by the phase of the LO, more precisely 

Z ( / >L (x,t) = a^{x,t)a b D {x,t) — a c ^{x,t)a c D {x,t) (14) 

(15) 
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where Pl(x) = |ax,(x)| and 4>l{%) = arga^x). Taking into account the finite size of the 
pixel detection area Rj, the measured quantity is 

zf L {t)=\ dxZ H {x,t). (16) 

JRj 

We now want to compare the fluctuations of the field quadrature measured in two symmet- 
rical pixels j = 1 and j = 2 of the signal and idler image. To this aim, we consider the sum 
and the difference of the quadrature obtained from two symmetrical detection regions R\ 
and R2: 

zf I !(t) = Z { l ) L {t)±zf L {t). (17) 
The corresponding fluctuation spectra, defined as 

V it\ Q ) = r dte™(5z£\t)5z£\0)) , (18) 
j— 00 

sz£\t) = Z%\t) - (Z%\t)) , 

describe the degree of correlation between the observables and Z^. Using the input- 
output transformation (Q), the commutation rules (0) and the fact that the input image is 
coherent we obtain the following relations: 

= (l-v)(SN) LO +v 2 [ dx\u(x,Q)al(x)-v*(x,Q)a L (-x)\ 2 , (19) 

JR1+R2 

where 

(SN) LO = vf dxp\{x) (20) 

is the shot noise level determined by the LO on the two detectors (we assumed |a^(x)| 2 ^> 
(j(x,t))). Next, we assume that the LO is symmetric with respect to the system axis, i.e. 
ocl(%) = oil{—x). Because v(x, Q) = v(x, —Q) we can write: 

Vj- L \Q) = V£l /2 (Q) = (1 - V )(SN) LO + v 2 I dxF(x, Q)p 2 L (x) (21) 

where 

F(x, Q) = \u{x, tye^^ - v*(x, fi)e^ w | 2 . (22) 

Note that the above expression corresponds to the fluctuation spectrum normalized to shot 
noise when 77 = 1 and the pixel area is small with respect to Xq and to the square of the scale 
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of variation of a(x). In this case x in Eq.(|2T|) must be taken as the central point of pixel 1 or 
pixel 2; the result is the same for both pixels because u(x, fl) = u(—x, fl),v(x, fl) = v(x, fl) 
and 4>l{x) = <f) L {-x). 

A first important result follows from the first equality (0)), according to which Z$ 

(2) 

and are correlated one to each other exactly to the same extent as the corresponding 
orthogonal quadrature components Z^ +7r , 2 and zjj^ +7r , 2 are anti-correlated. Second, the 
common fluctuation spectrum of the two observables Zjj^ — Z^ and Z^ +7r , 2 + Z^ +7V , 2 as 



given by expression (|T9| ) does not depend on the intensity and phase of the input image. 
Hence the result is the same in the phase-insensitive and in the phase-sensitive scheme, and 
remains the same even in absence of an input image at all, i.e. in the case of pure parametric 
fluorescence. Third, this spectrum can be reduced well below the shot noise level, provided 
the gain is large enough and the phase of the LO is correctly adjusted. Indeed, assuming 
that 

<t>L{x) = ^ (arg w(£, 0) + axgv(x, 0)) = (j) opt (x) (23) 

over the two detection areas, using the symmetry property of the LO oll{— x) = ot^x) and 
unitarity relations ([|), one obtains from Eq. (|T9|) for Q = 0: 

F(x, Q = 0) = ? , (24) 

[Kf,o)| + Kf,o)|] 2 

which goes to zero when \u(x, 0)| ~ \v(x, 0)| ^> 1. Under conditions of large gain and reason- 
ably large quantum efficiency, almost perfect correlation between the selected quadratures 
can therefore be obtained. 

It is interesting to relate the phase of optimum squeezing (in Z$ — Z$) 4> op t for the LO 
with the phase of maximum amplification in the phase sensitive configuration. The mean 
output field is in general 

<y ou t{x) = (a 4 {x, t)) = u(x, 0)a in (-x) + v(x, 0)a* n (x), (25) 

where we used Eq.(^) and Eq.(||). In the phase sensitive case c*m(— x) = a» n (x) we can write 

\a out (x)\ 2 = G(x)\a m (x)\ 2 (26) 

with the phase-sensitive gain given by: 

G{x) = \u(x, 0)e^^ + v(x, o) e -^"^| 2 , (27) 



where <fr% n (x) is the phase of di n {x). We easily obtain that the maximum gain 

G max {x) = \\u(x,0)\ + \v{x,0)\\ 2 (28) 

is obtained for: 

<f>in($) = 4>?n x {$) = \ (argU(f, 0) - argU(f, 0)) . (29) 



L xmax 

Hn — Yin 



On the other hand, from Eq.(^) with ««„(— x) = ai n (x) one obtains that when 
the phase 4> out of the output field a out is given by 

<Pout(x) = <F2£{x) = ~ (argu(f , 0) + axgv{x, 0)) (30) 

and therefore coincides with op t(x) given by Eq. (|23|) . This leads to the following interpre- 
tation for <j) opt : the phase 4>l of optimum squeezing in Z^T: coincides with the phase of the 
output field in the phase-sensitive configuration, provided the phase of the input field is se- 
lected to have maximal amplification. Note that in the special case of perfect phase matching 
one has &rgu(x, 0) = so that <fi™ x (x) given by Eq.(^) coincides with the corresponding 

4>out(x). 

The results obtained for the observables Z^r closely ressemble the situation of the EPR 
paradox for continous variables demonstrated in || , but generalised to many pixels (see also 
0) and to the presence of input images. We notice indeed that the conjugated observables 
Xj = f_^\dtZ { l\i) and Pj = jf^/a d *^??-Hr/ 2 (*) U = M) obe y the uncertainty rule: 



-T D /2 UA ^ J <f> L \. v J al1 ^ J - J-T.d/2 """"fo-Mr/S 

1 r i2 



T D / dxp 2 L (x) 



(31) 



IR1+R2 

On the other hand, the following combination over the two pixels: X_ = X% — X2 and 
P + = P\ + P2 are commuting observables that can be simultanously determined. When 
the time of measurement To is much larger than the inverse of the temporal bandwidth of 
the OPA, using Eq. (|I8D , the uncertainty of these observables can be directly related to the 
fluctuations spectrum 

(5 2 X^) = (5 2 P + ) = T D v£(n = 0) . (32) 

For 1] = 1, an optimal adjustement of the LO phase allows these uncertainties to reach 
almost a zero value for large amplification and thus to display an apparent violation of the 
Heisenberg rule: 

1 r i2 



(5 2 X„)(5 2 P + } < 4 



T D [ dxp\(x) . (33) 
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However, it is impractical to synthetyze a LO with the phase variation prescribed by 
Eg. (p3|) . On the other hand, for a LO with constant phase, the condition ( |23|) concerning 
the phase of the LO can be exactly satisfied only for a single couple of pixels of area small 
compared to So, so that the gain functions \u(x, 0)| and \v(x, 0)| are nearly uniform over the 
detection areas. We can however show that , by introducing an appropriate curvature the the 
wavefront of the LO field, EPR-like correlations are present for each couple of symmetric 
pixels in the output over the whole gain region Sq. To this end, we allow the LO phase 
distribution to have a quadratic dependence on the spatial coordinate (which corresponds 
to a spherical wavefront as one has e.g. in gaussian beams). The wavefront curvature is 
selected in order to have the best fit of the spatial dependance of (f> op t(x) in Eq . (|23|) . 

Figure 2 plots the function F(x,Q = 0) in the limit where R\ and R2 are small compared 
to Sq and symmetric. The collinear phase-mismatch at degeneracy is A / c = 0.5 and the 
linear gain parameter is \a\l c = 1.5. 

Curve (a) corresponds to the ideal case, with satisfying condition ( p3|) everywhere 
in the transverse plane, and leads to a maximal amount of noise reduction in the whole 
amplification region. In curve (b), the phase of the LO is constant and satisfies condition 
(p3|) only in the point of maximum gain xq, where perfect phase-matching is achieved, 
\u\ 2 ~ 5.5 and F = (\u\ — \v\) 2 = exp(— 2|<r|Z c ) ~ 0.1. Curve (c), obtained by optimizing the 
phase with a quadratic term (i.e. we take the form 4>l(x) = $ + ^(l^l — x g) 2 / x o ), is the 
best that can be done with a gaussian LO and is close to the ideal case. 

IV. PHASE-INTENSITY ENTANGLEMENT OF THE TWIN IMAGES 

Although the phase-sensitive measurement scheme considered in the last section offers 
a picture of the spatial correlations that can be observed in the output field, intensity 
correlation measurements are more straightforward to perform experimentally and lead also 
to interesting effects of quantum noise reduction [Q, f|, [14|. The observable that displays 
reduced fluctuations with respect to the coherent state level is the difference between the 
direct photocurrents measured from two symmetrical detection region z_ = i\ — i 2 . The 
corresponding fluctuation spectrum is 

/oo 
dt e int (5i^(t) <M0)> (34) 
-00 
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By using (Q), @, (||), (|TOl), ( pTTj ) and the fact that input image is in a coherent state one 
obtains after lengthy but elementary calculations 



VL(fi) = (l-*7)<i+> 

+?7 2 £ df |w(f, - if (£, -Vt)a out (-x)\ 2 

\u(x, n')\ 2 \v(x,n + n')\ 2 
-m(x, n + q')v(x, -n - n')u* (x , n')v* ($, -n')] (35) 



R1+R2 
2 1 r ^ r dQ! 



+77 — / 

JR1+R2 J 



where a 0?1 t(af) is given by Eq. (|25"D . The shot noise level corresponds to the photocurrent sum 
= {ii+i2}- The second term on the l.h.s. ofEq. (|35|) arises from the interference between 
the amplified input field and the fluorescence field. The last term, which does not depend 
on the presence of an input field, is a pure noise contribution due to the self-interference of 
the fluorescence field and reduces to zero for Q = because v(x, —Q) = v(x, £1) [0. 

Using the explicit expression of the amplified input field ( p5|) and the fact that u(x, Q) = 
u(—x, Q), v(x, fl) = v(—x, fl) and Eq.(H), we find for the zero frequency value of the spectrum 

V i _(0) = (l-r l )(i + )+r l 2 [ dx\a m (x)\ 2 (36) 

JR1+R2 

As shown in in the case of ideal detection (rj = 1), the noise level of i_ reduces therefore 
to the noise of the input image over Ri + R 2 . As a consequence, under conditions of large 
gain, fluctuations are well below the shot noise level. 

It is important now to connect with the result for quadrature components obtained in 
the previous section. To this aim, let us first assume that the input field is strictly different 
from zero at least in some region of the transverse plane. Second, let us assume that the 
parametric values are such that the pure noise contribution in (Q) (i.e. the last term in 
Eq.(pS)) is negligible and, similarly, that the second term in Eq.(|ll|) can be dropped. Thus 



expression fl3"oD) reduces to 

v^n) = (i- v )SN- 

+r] 2 dx \u(x, VL)a* out (x) -v*(x, -Vt)a out (-x)\ 2 , (37) 
JR1+R2 

where 

SN_ = r][ dx \a out {x)\ 2 (38) 

JRi+R?. 



IR1+R2 
12 



F(x, n) = \u(x, n)e- l ' t,0 ^ ) - v*{x, -fi) e ^°" t(:?) | 2 , (40) 



and we used Eq.(^). By comparing with Eqs. (|19"D and (0), we see that this expression 
coincides with V^ L if we take: 

a L (x) = a out {x) . (39) 

This is expected because a LO with the configuration of the output field just picks up the 
amplitude fluctuations. The relation becomes even more precise in the phase-sensitive case 
&in{x) = oci n (—x). In this case, assuming rj = 1 and that the pixel area is small with respect 
to Xq and to the square of the scale of variation of a ou t(x), one has: 

SN_ 

where we set a out (x) = p out (x) exp(icj) out (x)), and x is the central point of any of the two 
symmetrical pixels. The result coincides with that of Eq.([H|) where 0l is replaced by (pout- 
The link between intensity fluctuations and quadrature fluctuations allows now to analyse 
immediately the case of phase fluctuations, which coincide with the quadrature fluctuations 
obtained by using a LO which displays a phase shift of it/ 2 with respect to LO which 
provides the amplitude fluctuations. Therefore, in terms of pixels, we are lead consider the 
observables Z^ +7T ^ 2 (t) (see Eq.(|T6|)) with <p L = 4> out . This naturally induces us to focus on 
the observable Z<+>(t) = ^I +7r/2 (*) + ^ O it+7r/2(0> see Eq.(PD) which measures the degree 
of anticorrelation between the phase fluctuations in the two symmetrical pixels 1 and 2. The 
spectrum V^^, 2 (p,) coincides with vj~\&), which as we have seen is identical to V^_(0) 
given by Eq.(^). Therefore for large amplification the fluctuations of Z^ + \t) are well below 
the shot noise level, which implies that the phase fluctuations in the two symmetrical pixels 
are strongly anticorrelated, exactly as the amplitude fluctuations are strongly correlated. 



V. CONCLUSION 

In this article we analyzed extensively a system formed by an optical parametric amplifier 
with some imaging lenses. Amplification of optical images by OPA has been already studied 



in the literature |T5[], but only from a classical viewpoint. 

Our results hold both for a phase-sensitive configuration (symmetrical input image) and 
for a phase insensitive one (asymmetrical injection). 

We demonstrated that the two output twin images exhibit a complete spatial EPR en- 
tanglement. This was shown, first of all, by considering a pair of orthogonal quadrature 
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components of the output field. In the case of local oscillator symmetrical with respect to 
the system axis, we found a precise prescription for the phase 4>l of the local oscillator (see 
introduction) in order to observe maximal correlation between symmetrical pixels of the two 
output images. The optimal value for the phase is that which corresponds to the amplitude 
fluctuations of the output images in the phase sensitive configuration, when the phase of the 
symmetrical input images is selected to obtain maximal amplification. 

We have shown also that a performance very close to that of the ideal case of the optimal 
LO phase can be obtained in practice by using a LO with a quadrature wavefront (as one 
has in gaussian beams) with the curvature used as optimisation parameter. 

The connection between quadrature fluctuations and amplitude/phase fluctuations has 
in turn allowed us to conclude also that, while intensity fluctuations are strongly correlated 
in the twin images, phase fluctuations are strongly anticorrelated in the same amount. An 
amusing analogy with amplitude and phase fluctuations in entangled twin images is provided 
by a fossile broken in two pieces (see Fig. 3). We see that the structures in the two pieces 
have the same "amplitude" = thickness, but opposite "phase" (one is concave and the other 
convex, one is righthanded and the other lefthanded). 

It is important to underline that, while the results for intensity and phase fluctuations 
hold only in presence of an input image, the result on EPR entanglement of quadrature 
components hold also in absence of any input image, i.e. in the case of the pure parametric 



down-conversion as in ||14|| . This is important for the applications to quantum teleportation 



of optical images ||16|| , as a generalisation of the Braunstein-Kimble [|17], |18| scheme for a 
single mode field, or to quantum cryptography with images. 

We observe finally that our results hold also if the OPA is replaced by an optical para- 
metric oscillator below threshold with plane mirrors (see || in this connection). As a matter 
of fact, also in this case one has an input-output relation of the form (|2|), and the results 
are based only on this relation and on the general properties of the functions u and v. 
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VI. FIGURE CAPTIONS 
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FIG. 1: Schematic diagram of the parametric image amplifier. The two- lens telescopic system 
allows to obtain two amplified copies of the input image that are strongly quantum correlated to 
each other, thereby the name twin images. The device is phase-sensitive when the input image 
is symmetrical, phase-insensitive when it is confined in the upper half of plane Pi. / is the focal 
distance of the lenses. 



FIG. 2: Plot of the noise reduction factor F(x, 0). Subscripts (a) refers to the optimal phase of 
the LO while (b) and (c) refer respectively to a constant phase and to a phase with a quadratic 
dependence on the distance from the optical axis. The dashed line is the phase-sensitive gain of 
the OPA (see Eq.@ divided by a factor 10). A l c = 0.5 and \a\l c = 1.5. 



FIG. 3: Analogy between a broken fossile and quantum entangled images (see text). 
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